SYNCHRONISATION PROPERTIES OF TREES IN THE 
KURAMOTO MODEL 
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Abstract. We consider the Kuramoto model of coupled oscillators, specifically the case of tree 
networks, for which we prove a simple closed-form expression for the critical coupling. For several 
classes of tree, and for both uniform and Gaussian vertex frequency distributions, we provide tight 
closed form bounds and empirical expressions for the expected value of the critical coupling. We also 
provide several bounds on the expected value of the critical coupling for all trees. Finally, we show 
that for a given set of vertex frequencies, there is a rearrangement of oscillator frequencies for which 
the critical coupling is bounded by the spread of frequencies. 
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1. Introduction. The Kuramoto model [161 [151 [71 [10] "was originally motivated 
by the phenomenon of collective synchronisation whereby a system of coupled oscillat- 
ing vertices (nodes) will sometimes lock on to a common frequency despite differences 
in the natural frequencies of the individual vertices. Biological examples include oscil- 
lations of the heart [3T] and of chemical systems [12]. Systems of coupled oscillators 
can also be used as an abstract model for synchronisation in organisations [31 IH H] . 
While Kuramoto studied the infinite complete network it is natural to consider finite 
networks of any given topology. This would correspond to a notion of coupling that 
is not universal across all vertex (node) pairs, but rather applies to a subset of all 
possible edges (links). For example the work patterns of human individuals in an 
organisation might enjoy a coupling effect in relation to pairs of individuals that have 
a working relationship. 

A typical network of coupled oscillators is as shown in Figure [U Each vertex has 
an associated phase angle 9i, as well as its own natural frequency w^. In this paper 
we are primarily concerned with the case where the network is a tree. The basic 
governing equation is the differential equation: 

n 

(1) ei=uji + kY^ A,j sin(6ij - Oi) 

1=1 

where A is the adjacency matrix of the network and A: is a coupling constant which 
determines the strength of the coupling. We refer to a network (graph) with preferred 
oscillator frequencies uji attached to the vertices as a Kuramoto graph. 

It has been observed [3 117) that many Kuramoto graphs synchronise^ in that the 
actual vertex frequencies 9i converge to a common value. That is, the vertex phases 
9i rotate at the same rate, with a constant phase difference between each pair of 
vertices. Moreover this phenomenon appears at a critical coupling constant kc, and 
as k increases, this phenomenon applies to a greater range of initial phases (the 9i at 
t — 0). Thus the graph has a frequency fixed point characterised by all the 9i being 
equal. 
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Fig. 1. An example Kuramoto tree. 

At the frequency fixed point, tlie Oi are equal to ui, tfie mean of the frequencies 
Wi, and it is convenient to apply a rotating frame of reference, with: 

(2) Mt) ^ 0,{t) - Cot 

We then obtain a system of differential equations equivalent to the original: 



(3) 



^ uji — uj -\- k Ai 

i=l 



At the frequency fixed point, we then have 0^ = for all i. While the literature 
contains several definitions of "critical coupling," here we define the critical coupling 
fee as that number for which a frequency fixed point exists exactly when k > kc (in 
[5], this critical coupling is called Kl)- 

In general the value of the critical coupling is the solution of simultaneous tran- 
scendental equations and would not be expected to have a closed form solution. How- 
ever, for complete graphs and complete bipartite graphs the critical coupling can be 
computed efficiently as the solution of non-linear equations [T51 [12] . For the graph 
with only two vertices i^i and 1^2 and a single edge between them, it is easy to see that 
solutions exist precisely when \uJi — < 2k, i.e. 



(4) 



kr 



IcJi - UJ2\ 



Our first theorem generalises this to any tree. We then use that result to calculate 
the expected value of the critical coupling for several classes of Kuramoto tree having 
randomly chosen frequencies, and to find bounds on the expected value of the critical 
coupling for Kuramoto trees in general. Finally, we use the result to show that 
the critical coupling for a specific tree topology can be reduced to a value which is 
independent of the number of vertices, by reshuffiing the frequencies w^. 

2. Main Lemma and Corollaries. The first step to our initial theorem is a 
lemma which simplifies the calculation of critical couplings by breaking a Kuramoto 
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Fig. 2. Construction for Lemma^^ 



graph into components. In particular, we consider Kuramoto graphs with a cut- vertex, 
i.e. one whose removal disconnects the graph. Figure [2] shows an example. 

Lemma 1. Let G he a Kuramoto graph with n vertices ui,i'2, ■ ■ ■ ,Vn with nat- 
ural frequencies Wi , a;2 , . . . , in which vertex Vi is a cut-vertex ( whose removal 
disconnects the graph) with vertices 1/2, ■ ■ ■ ,Vm on one side of the cut and vertices 
Vm+i, ■■■T^n on the other. As shown in Figure\^ let Gi and G2 be the two graphs 
formed by deleting all the vertices on one side of the cut vertex and replacing the cut 
vertex by new vertices x (for Gi ) and y (for G2 ) having frequencies lUx and ujy : 



(5) 



UJ ■ 



i=2 



(uji — Uj) — UJl + 



i—m+l 



UJ 



(6) ujy = UJ — ^ (ujj — uj) = UJl + y^(a^i — uj) 

i—va+l i—2 

Then G has a frequency fixed point if and only if both Gi and G2 have frequency fixed 
points, and hence kciG) — max{kc{Gi),kc{G2))- 

Proof. First we note that the frequencies uJx and ujy are chosen in [5] and [S] so that 
the average frequencies of both Gi and G2 are the same as that of G. Let G have a 
frequency fixed point so that we have k and (f>i,i = 1, . . . , n satisfying the system of 
equations for i = 1, . . . ,n: 

n 

(7) = Wi - w + fc ^ Aij sm{(pj - 

i=i 

If we sum these equations from i — l,...,m, then each term Aij sm{(j)j — (pi) is 
cancelled by the term AjiSm.{(j)i — with the exception of the pairs = 
m + 1, . . . , n and we obtain: 



(8) 



m n 

= ^(wi-a;) + fc ^ Ausm{4>i - 

i—1 -i—m+l 
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Subtracting equation ([S]) from ([7]), for the case i = 1, gives an equation for vertex x: 



i=2 



(9) 



= - ^(wi - w) + fc ^ Ah sin(0j 

! i=2 
m 

Ox + An sin((?!)i - 0i ) 



i=2 



Similarly by summing ([7]) from i = m + 1, . . . , n we obtain an equation for vertex y: 



(10) 



= cj — oj,, 



k y: a 



li sm( 



i) 



It follows that a solution to the system of equations ([7]) is a solution to the system of 
equations for Gi and G2 in which the vertices x of Gi and y of G2 have the same phase 
(j)i . If on the other hand we have a solution to the system of fixed point equations for 
Gi and G2, then we can shift all the phase angles of the vertices of G2 by the same 
amount so that (j)y — without changing the phase angle differences. This provides 
a fixed point solution to the system of equations ([7]) and completes the proof. □ 

Repeated application of this lemma can reduce the calculation of the critical 
coupling to 2-connected graphs (in which any two vertices lie on a cycle) and tree 
edges. For example, we have the following corollary for graphs with a cut-edge, i.e. 
an edge whose removal disconnects the graph: 

Corollary 2. Let G be a Kuramoto graph with n vertices 1,2, ... ,n with natural 
frequencies loi,uj2, . . . ,ujn in which the edge from i>i to i>2 is a cut-edge with vertices vi 
and 1^3, ... , Vm on one side of the cut and vertices V2 o-nd Vm+i, ■ ■ ■ i^n on the other. 
Let G' be the one-edge graph containing only the vertices x and y, where: 



(11) 



(12) 



1=3 



UJy = UJ2+ / 

z— rn+l 



Then G has a frequency fixed point only if G' has a frequency fixed point, i.e. 



(13) 



^(Wi -Cj) - 



E 



(Wj - Uj) 



< 2k 



Proof. Applying Lemma [T] twice yields three components: one based on the 
vertices {vi,v^, . . . ,Vm\, one based on {v2,Vm+i, ■ ■ ■ ,Vn\, and the one-edge graph 
G' = {x,y} based on {vi,h'2}. The condition on k follows by applying (U) to G'. 
□ 

For graphs which are trees, repeated application of Lemma[l]produces components 
which are exactly the tree edges, with a similar adjustment of frequencies. We can 
therefore obtain an "if and only if" condition, and hence obtain a precise value for kc- 

Theorem 3. Let T be any tree with at least two vertices. For any edge e ofT, let 
Tf, and be the two subtrees formed when e is deleted from T . These two partitions 
have |Te| and \T'J^ vertices respectively. Define the partition sum: 



(14) 



0(T,) = 
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be the sum of the frequency deviations from oj. Then ^1{Tq) — ^{T^) and the critical 
coupling for T is given by: 

(15) kc = maxr2(Te) 

Proof. We repeatedly apply Lemma [1] to obtain n — 1 graphs each consisting of 
a single edge. The condition on k follows from (jl)), applied to each of the n — 1 edge 
graphs obtained. □ 

We note that, for trees, the condition in this theorem is equivalent to the more 
general sufficient conditions given in Remark 10 of [8] and Statement Gl of [6]. In 
the remainder of the paper we consider various implications of the result, facilitated 
by the form in which our result is expressed. We begin by considering the expected 
value of fee, when frequencies are chosen randomly. 

3. Expected Critical Couplings. It may be the case that we do not know 
the precise frequencies for a Kuramoto tree, but we do know the distribution of the 
frequencies. In this case, we can apply Theorem [3] to determine the expected value of 
the critical coupling kc for various kinds of tree. Figure [3] shows some of the trees 
considered. In this section we combine theoretical derivations, using Theorem [3] and 
statistical theory, with Monte Carlo experiments. In these experiments, we consid- 
ered a variety of trees, and made 1,000,000 frequency assignments with uji uniformly 
distributed over the interval [0,1]. The frequencies oji therefore have mean 1/2 and 
variance ct^ = 1/12, and the expected critical coupling E{kc) is proportional to cr, 
the square root of the variance. We used Theorem [3] to determine the critical cou- 
pling kc for each assignment, and hence the expected value of the critical coupling for 
that size and shape of tree. We also performed the same experiments using normally 
distributed frequencies with the same mean and variance, calculating over 8,000,000 
frequency assignments in that case, to allow adequate sampling of the tail. 

3.1. Chains. For the case of a chain (a tree with all vertices connected in a line) 
with n vertices vi, . . . ,Vm by Theorem [31 the critical coupling kf^^^^ is: 



(16) fcf = max 



^(w, - w) 



The partition sums X]i=i(^i ~ ^) f*^^ ^T:st j vertices form a random walk as j 
increases. Figure |4] illustrates an example, for n = 60. The walk must return to 
zero, since u} — Y^^=i ^^'^ the maximum translation distance is likely to occur 
in the vicinity of the centre of the chain. This random walk has step variance cr^, 
and a crude estimate for the expected critical coupling E{kf^'^^'^) is the expected final 
displacement A for a random walk with this step variance and r steps, where r is 
some value between n/2 and n. However, this takes into account neither the return 
to zero nor the fact that the maximum in (jl6l) is being taken over all values of j. An 
asymptotic expression for A was given by CofFman et al. [3 [5] which, for the case 
T = n/2, gives: 




(17) A = (T\l^^ - (yc = a\j'- -ac = -0.148976 

where c = 0.516068, and we have ignored terms in n~^/^. Although this expected 
translation distance under-estimates the critical coupling, we will still have a result 
of the form a^/n + 6, and we can determine the values of a and h empirically. 
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Fig. 4. One random walk of partition sums X]i=i('^i ~ chain with n = 60 and 

increasing j . Dotted lines show the empirical expected value E{k'^^^^'^) = xC"-); which is the expected 
peak deviation from 0. 



For our empirical Monte Carlo studies, we considered chains with the number 
of vertices n ranging from 2 to 300, with 1,000,000 different uniformly distributed 
frequency assignments (and 8,000,000 different normally distributed frequency assign- 
ments) on each. We used Theorem |3] to calculate critical couplings in each case, and 
fitted a curve of the form ay^ + b to the expected values. Empirically, as illustrated 
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in Figure [SI the expected critical coupling for a chain of n vertices is: 

(18) £;(fcf ^'") = xin) ~ 0.252V^ - 0.168 = 0.873 aV^ - 0.581 a 

We express the result both in terms of a, and for our special case of cr = y/l/12. As 
we will see in WSM the expected critical coupling for the chain acts as an upper bound 
on E{kc) for trees in general. To facilitate use in this way, we abbreviate this function 
as x(ri). An upper bound on x('i) itself is the expected maximum displacement for n 
steps given by Weiss [20l p. 192]. This is too high because it ignores the constraint 
of returning to zero: 

(19) E{kt^'^)^x{n)<a^^ 

Figure [5] shows this bound, together with the corresponding lower bound, which is 
the expected maximum displacement for n/2 steps. This is too low because it ignores 
the fact the peak deviation from zero can occur at either end of the chain: 

(20) E{kf-'^^)^x{n)>a^^ 






Order n 



Fig. 5. Expected critical couplings for the chain, as a function of the number of vertices n. 
The solid line shows the empirical curve 

^(•^jchain-) ^ ^ 0.252^^-0.168, while the dashed and 
dotted lines shows the random-walk bounds \J ttu / 48 and ^JwnJ^. Results for normally distributed 
frequencies are, as expected, virtually identical, and cannot be distinguished on this plot. 



3.2. Stars. For star trees, such as the one in Figure[3](a), where i^i is the central 
vertex: 

(21) /cf" = max |wi -w| 

i—2...n 

For our choice of frequencies, taken from a distribution with mean 1/2, the critical 
coupling, which is the expected value of this maximum, is given approximately by: 

(22) Eikf""') w E fmax 
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where the second term is small. The expected maximum of — which is uni- 
formly distributed, is given approximately by: 



(23) 



E ( max 



1 

UJi 

2 



n - 2 



n- 2 
2n 



This uses the fact that the expected maximum of n independent uniform distributions 

is m- 



n- 1 
n+1 



and the fact that the distributions are partially correlated, so that one of the uJi can 
be inferred from uj and the other Ui. In addition, the expected value of ui differs from 
1/2 by approximately: 



(24) 



E 



1 



1 

6mr 



By ([22l), ([231), and a good estimate for the expected critical coupling E{kf'^'') 

will be: 



(25) 



n-2 



2n 



1 

6n7r 



Figure [5] shows that, for n > 40, this is an excellent estimate of the actual critical 
couplings, which are shown by solid triangles. It can also be seen that, as n tends 
to infinity, the value of E{kf^'^) approaches 1/2. The expected critical coupling for 
the star acts as a lower bound for trees in general, since star trees have the smallest 
possible partitions. 

For comparison, with normally distributed frequencies, the expected maximum of 
[wi — ai| can be very closely approximated using the expected value of the maximum 
of the absolute values of x independent normally distributed random variables of unit 
variance, which we write as /i(a;). For a; = 1, /i(l) = ^21'k. For a; > 2, the value of 
Pu{x) is slightly higher than (and converges to) the expected value of the maximum of 
X independent normally distributed random variables of unit variance and zero mean, 
which can be roughly approximated by •\/2 log(a;) |llj . The function [i(x) can be 
expressed more precisely using the inverse complementary error function: 



(26) 



\ji{x) = 72 



(1-7) erfc"^ 



- ) + 7 erfc ^ f — 
X \xe 



where 7 — 0.57721566490 ... is the Euler-Mascheroni Constant. However, we set 
/i(0) = by definition. 

The expected critical coupling for normally distributed frequencies, shown by 
open triangles in Figure [51 is very closely approximated by: 



(T/i(n — 2) 



fi{n - 2) 
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Fig. 6. Expected critical couplings for the star, as a function of the number of vertices n. The 
thicker solid line shows the curve (n — 2)/ (2n) + ^l/(6ra7r), while the thinner line shows the infinite- 
size limit of 1/2 for uniform distributions. For comparison, open triangles show expected critical 
couplings for normally distributed frequencies, and the dashed line shows the estimator crfi{n — 2). 




Number of standard normals x 



Fig. 7. The function tJ,{x), the expected value of the maximum of the absolute values of x 
independent normally distributed random variables of unit variance (solid line). The dashed line 
shows the approximation ^2 log(x), while the dotted line shows fJ.{l) = \/2/tt. 



3.3. Dumb-bells. For "dumb-bell" trees, such as the one in Figure |3] (b) , the 
critical coupling is the result of the combined effects of the n — 2 leaves and the fact 
that the tree can be partitioned into two halves. The latter has the greater influence, 
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SO by Theorem [31 
(27) 



n/2 




n 


^(w, - w) 






i=l 




i=l+n/2 



That is, the critical coupling fc^ is the deviation from uj for the average of the uji for 



one half of the dumb-bell (without loss of generality, we will use uji , . 
calculate E{kf^), the expected value of this deviation, we note that: 

(28) 

Thus: 
(29) 



i/2 



). To 




n/2 

E 



n _ 

— — LlI 



n/2 ^ n ^ n/2 

E^»~ 2^^* ^2^"^' 

i—1 i—1 i—1 



1 " 

i=l+n/2 



The variance for this difference of sums is and hence the expected critical 

coupling will be: 



(30) 



m 



n 
24^ 



The dashed line in Figure [8] shows this approximation. A better approximation can 
be obtained by noting that when calculating the critical coupling from |wi — wj for 
half the dumb-bell gives a result of less than 1/2 = o-s/S, then the critical coupling 
will be determined by the effect of the leaves, in a similar way to the expected value 
for the star. Since the variance for (|29p is 7ict^/4, the standard deviation is: 



(31) 



Using a first-order approximation, the probability that a normal distribution with 
that standard deviation is within 1/2 = a\/2i of the mean will be approximately: 

o cr\/3 



crV27r 

The impact of the leaves is, with this probability, to replace an expected value of 
about 1/4 = cr(\/3)/2 by one of about 1/2 = <t\/3, i.e. to increase the expected 
critical coupling to: 



(32) 



Figure |8] shows that this correction ensures an excellent fit for n > 20. Note that as 
n increases, the correction factor tends to zero, and E{k'^^) tends to the value given 
by dSO])- 




SYNCHRONISATION PROPERTIES OF TREES 



11 




o _ 
o 

^ 1 1 1 1 \ 1^ 

50 100 150 200 250 300 

Order n 

Fig. 8. Expected critical couplings for the dumb-bell. The dashed line shows the initial approx- 
imation ^n/(247r) and the solid line the improved approximation, which adds the correction factor 
^3/(2n7r) due to the leaves. For comparison, open squares show values for normally distributed 
frequencies, which are higher due to an increased contribution from the leaves. 

3.4. Binary trees. We also considered binary trees, such as the one in Figure 
131(c), i.e. trees with a root vertex vi, and with vertices i>2i and 1^21+1 having Vi as a 
parent vertex. For these trees, the calculated expected critical couplings were well- 
predicted by the empirical curve 0.212-y/n — 0.082, shown in Figure [9l A precise 
analytic solution for this case is difficult, since the binary tree contains partitions of 
multiple sizes, ranging from 1 to (n — l)/2. However, as Figure |9] shows, the binary 
tree is intermediate between the chain and the dumb-bell. 

3.5. Tadpoles. "Tadpole" graphs, such as the one in Figure 12] (d) , consist of a 
star combined with a "tail," so that the diameter D > 2. Essentially, this is a kind of 
asymmetrical dumb-bell. The results in Figure [TOl are based on tadpoles with D = 8. 
For our uniformly distributed frequencies, the critical couplings converge to a limit of 
about 0.85, which reflects the combined contributions of the tail and the star. Such 
a combination of effects is typical of non-symmetrical trees. For normally distributed 
frequencies, the expected critical coupling grows logarithmically, as with the star. 

3.6. Upper and lower bounds. We now consider bounds on the expected 
critical couplings for trees with uniformly distributed frequencies. It is helpful to 
define as the number of vertices in the smaller of the two sub-trees Tg and 
formed when e is deleted from T, and to define P as the maximum of the Pg- We call 
P the maximum partition .size. 

In terms of n, the chain provides an upper bound for expected critical couplings 
on trees, since rearranging the topology of a chain can only result in some or all 
of the Pe being smaller. In terms of diameter, however, the chain provides a lower 
bound, as Figure [11] illustrates. This follows from considering the longest subchain 
(with D -\- 1 vertices) within a tree. Any additional vertices can only increase the 
expected critical coupling. Consequently, the critical coupling for a tree of diameter 
D is bounded below by the critical coupling for a chain of diameter D. Since x{^) 
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Fig. 9. Expected critical couplings for the binary tree. The solid line shows the empirical curve 
0.212-\/ri — 0.082, For comparison, the dashed line shows the empirical curve for the chain, while the 
dotted line shows the theoretical curve for the dumb-bell. Results for normally distributed frequencies 
are, as expected, virtually identical, and cannot be distinguished on this plot. 
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Fig. 10. Expected critical couplings for the D = 8 tadpole, as a function of the number of 
vertices n. For our uniformly distributed frequencies, the critical couplings converge to a limit of 
about 0.85. For comparison, open circles show expected critical couplings for normally distributed 
frequencies. These grow logarithmically, as with the star. 



gives the expected critical coupling for a chain of x vertices, we have: 
(33) xiD + l)<E{kc)<xin) 



For example, for a tadpole tree with D = 8, n > 17, and our choice of frequencies. 
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X(9) = 0.59 < 0.85 < 0.87 = x(17) 
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Fig. 11. Expected critical couplings for trees in terms of the diameter D. Downward-pointing 
black triangles represent trees not discussed above (including Y-shaped trees, X-shaped trees, and 
Scale-Free trees). The solid line gives the empirical curve x{-D + 1) for the chain. 

The maximum partition size also produces a lower bound. While the chain and 
the dumb-bell have different diameters, they both have P = [n/2\. For the binary 
tree, P = ['T./2J, while for the star, P — 1. For tadpole trees, P = D — 1, provided 
n > 2{D - 1). 

Partitions of size at least P must occur at least twice, and these partitions are 
not independent, since the uji sum to nui. Suppose that we decompose the tree into 
r partitions (each of the form for some edge e), that Xj is the sum of the cji over 
the j*"^ partition, and that the j"^ partition contains mj vertices. Then: 



(34) 



1, 



(1 - -)X 



Since we are considering partitions of size at least P, by a standard result on the 
variance of linear forms in distributions 13 , Xj 
at least: 



rrijUJ has a standard deviation of 



(35) 



a" = ay/P\ 



r - 1 



+ {r-l)\=<7y/pJl-- 



Since the number of partitions r > 2, we obtain a lower bound on E(kc): 
(36) 



^ /2 /P P P 

E{K) > <j\ -\h: ><^\ - = \hr 

V TT V 2 V TT V 127r 

Figure [T^ illustrates this lower bound, using the same datapoints as Figure [TT] For 
trees with small diameter, like the dumb-bell, this is tighter than the bound based on 
D. 
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Max partition P 



Fig. 12. Expected critical couplings for trees in terms of the maximum partition size P. The 
solid line gives the lower bound \J P/ (127r) . Datapoints are the same as Figure FTTI 



We also obtain an upper bound in terms of P by noting that the worst-case 
partition is half of a chain of length 2P, and that this partition occurs at most n/P 
times. We can approximate ^{n/P) by y^2\og{n/ P), so that, for suitable constants 
a, b, and c: 



/ 7^ — TL Tl 

(37) E{k,) < ax(2P)^2 log - = baVP^2\og - - ca^Plog - 

Since the chain is the worst case, from ([T8)) it suffices to take c = | > 0.873-\/2/ log 2, 
i.e. 



(38) Eik,) < ^a^Plogp = ^ySPlog^ 

Figure [T^ illustrates this bound. 

3.7. Bounds for normally distributed frequencies. Since normally dis- 
tributed frequencies lead to critical couplings that are either virtually identical to 
or larger than the critical couplings for uniformly distributed frequencies, and since 
the chain continues to be the worst case, the bounds in the previous section apply 
also for normally distributed frequencies. In particular. Figure [14] confirms visually 
that the upper bound ([38]) still holds. 

In terms of the order n, expected critical couplings are bounded above by x(^)) 
the value for the chain, as is the case for the uniform distribution. Expected criti- 
cal couplings are bounded below by afi{n — 2), the value for the star, which grows 
logarithmically. 

4. Minimising the Critical Coupling using Rearrangement. The results 
above have application to organisational trees. In this case, the frequencies uji can 
reflect the speed of the decision cycle for person i (as in 9 ) or, more abstractly, 
they can represent a one-dimensional projection of the activities of person i (as in 



SYNCHRONISATION PROPERTIES OF TREES 



15 





Fig. 14. Expected critical couplings for trees with normally distributed frequencies in terms of 
Plog(n/P). The solid line gives the upper bound (l/4)\/3Plog(n/P), as in Figure [TSl 



[5]). In both cases, synchronisation problems arise when vertices with different Wi 
are not strongly coupled. Organisational structures are often hierarchical trees (e.g. 
binary trees), because such trees offer a number of advantages, such as low diameter 
(worst-case distance between pairs of vertices) together with low degree (the number 
of subordinates an individual has responsibility for) [14] . If we also consider synchro- 
nisation to be favourable (in that steady organisational work patterns are achieved) , 
we would like to combine an efficient tree structure with a low critical coupling. 

As an example, Figure [K] shows the structure of a real organisation, where the 
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numbers in each vertex are a one-dimensional projection of each person's activities, 
calculated by administering a survey and taking the most informative of the principal 
components of the answers. Interpreting these numbers as frequencies, using Theorem 
[3] we calculate a critical coupling — 5.08. Improving the synchronisation of the 
network of people within this organisation requires either management activities to 
more strongly couple people's work, or some form of structural reorganisation. 




Fig. 15. An organisational tree with a critical coupling kc = 5.08, which reduces to 1.07 after 
rearranging frequencies. The two branches at the lower left are somewhat out of step with the rest 
of the organisation. 

As an example of structural reorganisation, for the organisation in Figure 1151 we 
can construct a Kuramoto tree which synchronises much more easily (with kc = 1.07) 
by leaving the topology of the organisational tree unchanged, but re-shuffling people 
(i.e. frequencies Wi) within it. In a managed situation like an organisation, we often 
have some discretion to reshuffle people in this way. In fact, the following theorem 
shows that we can reduce the critical coupling to a value which is independent of the 
number of vertices: 

Theorem 4. Let T be any tree network of n > 2 vertices, and uJi, for i — 1, . . . , n, 
be the collection of natural frequencies. Let T" be the tree which is isomorphic to T in 
which the frequencies are assigned to vertices according to some assignment a. Then 
for some a, 

(39) A:e(r°) < 

^max ^min 

where Wmax is the largest frequency and uJmin is the smallest. 

Proof. The proof is constructive. First we note that by traversing the edges and 
vertices of the tree in depth-first order, the corresponding assignment of the vertices 
of T to the natural numbers 1, 2, . . . , n has the property that for any edge e, at least 
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one of the components of T — e is assigned a set of numbers that are of the contiguous 
form fc, fc + 1, . . . , m for some k and m. In Figure [T51 we show an example of such a 
vertex assign ' ' ' ' ^ ' ' ' ' ' ous number 

sets. 




e,:{l}, e,:{l,2}, ei.{4}, e,:{l,2,3,4}, e,:{6}, e,:{7,8}, e,:{8} 

Fig. 16. An assignment of the numbers 1, . . . , 8 to vertices of a tree such that, for any edge e, 
at least one of the components ofT — e is assigned a contiguous set of numbers, as per Theorem \4\ 

We shaU assume that we have such an assignment of numbers to vertices and 
that vertex Vi refers to the vertex assigned to the number i. We now order the 
natural frequencies as follows. Select the first frequency to be Wmax- Add frequencies 
WQi,a;a2, • ■ • T^ap from any previously unselected frequencies less than uj, stopping at 
the least p where: 



(40) 



< 



pu) 



Add frequencies i^a{p+i)i'^a(p+2)j ■ ■ ■ i^aq from any previously unselected frequencies 
greater than uj stopping at the least q where: 



(41) 



9 



4=1 



LOai > 



In similar fashion continue to add vertices alternating between runs of frequencies less 
than w and greater than ui, until all frequencies are used. The function f{k) given by: 



(42) 



fik) 



.i=l 



must then take values between Wmax — ^ and Wmin ~ ^, with f{n) = (see Figure [TT]). 

Finally we assign frequency ujai to vertex Ui for i — 1, . . . , n. Recall that by our 
earlier observation, for any edge e at least one of the components of T — {e}, say Tg, 
has vertices with natural frequencies from a contiguous sequence uj„i, ^m+i, ■ ■ ■ , ^m+r 
for some m and r. These frequencies correspond to a segment of the function / from 
/(to) to /(to + r). Clearly: 



(43) 



|/(m + r)-/(m)| < (c^„ 



t^) - {uJmin -Uj) =UJn 



Also: 

(44) \f^m + r)-f{m)\^ 



m+r 



y2 ^ai - (r + \)u 



UJ.I - |Te|w 
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Combining (|43| and (j44|) . 



(45) 



Since this is true for all edges e, by Theorem [3] we obtain the result. □ 

The following example shows that this result is in a certain sense best possible. 
Let a star tree have 2 vertices with frequency ^ and n — 2 vertices with frequency C. 
Then for every assignment a of frequencies to vertices, 



(46) 



2C-2e 



le-ci 



This follows from Theorem|3l noting that every leaf vertex of the star must be assigned 
frequency ^ or 

Figure [18] shows an example for a 15-vertex binary tree. As another example, 
for 255- vertex binary trees, the expected critical coupling without rearrangement was 
3.289, using 1,000,000 different frequency assignments, with all frequencies uniformly 
distributed over the interval [0,1] (i.e. Wmax — Wmin < !)• Rearrangement ensured that 
all 1,000,000 critical couplings were restricted to the range 0.462 to 0.960 (less than 
"^max — "^min), with an cxpected critical coupling of 0.685. 

5. Discussion. In this paper we have explored the critical coupling for Ku- 
ramoto oscillators arranged in tree topologies. With Theorem [31 we provided a 
closed- form solution (|T5]) for the critical coupling fee, in the case that frequencies 
were known. 

In the case that only frequency distributions were known, we calculated Eikc), 
the expected value of the critical coupling, for chains and stars (for both uniform and 
Gaussian vertex frequency distributions) and for "dumb-bells" (for uniform distribu- 
tions only). We also provided empirical values of E{kc) for these cases, and for the 
"tadpole" and binary tree. For trees in general, we provided lower bounds on E{kc) 
in terms of the diameter D (|33p and the maximum partition size P (|36p . We provided 
upper bounds in terms of the number of vertices (the expected critical coupling for 
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Fig. 18. A binary tree with frequencies uniformly distributed over the interval [0,1], after 
rearrangement. Vertices are labelled with depth- first- search index and frequency The mean 

frequency ui is 0.527, and the critical coupling of 0.596 is based on the partition on the lower right. 

the chain performs this function) and the maximum partition size psp . These bounds 
hold for both uniform and Gaussian frequency distributions. 

FinaUy, Theorem 2] showed that, for a given set of vertex frequencies, there is a 
rearrangement of osciUator frequencies for which the critical coupling is bounded by 
the spread of frequencies. 

In future work, we hope to provide closed-form expressions for the expected crit- 
ical coupling E(kc) for a number of other cases of tree and, more importantly, to use 
methods similar to Lemma [T] to determine the critical coupling kc for other classes of 
graph. Corollary [5] is a small step in this direction. 
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